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Abstract 


The  goal  of  this  project  was  to  develop  new  and  innovative  processing  methods  for  clutter  and  interference 
mitigation  for  Space-Time  Adaptive  Processing  (STAP)  with  circular  arrays.  The  focus  was  on  developing 
a  robust  minimum  mean  square  error  (MMSE)  beamforming  technique  using  quadratic  quiescent  pattern 
constraints  (MMSE-QPC)  that  works  with  arbitrary  array  configurations,  including  the  circular  UHF  Elec¬ 
tronically  Scanned  Array  (UESA)  and  the  standard  linear  array.  The  MMSE-QPC  technique  provides  a 
general  framework  for  quiescent  and  adaptive  space-time  beampattem  synthesis  which  provides  both  main 
beam  and  sidelobe  control  with  reasonable  computational  complexity.  Beampattem  control  is  achieved  by 
imposing  a  set  of  inequality  constraints  on  the  weighted  mean-square  error  between  the  adaptive  pattern 
and  a  desired  beampattem  over  a  set  of  angle-Doppler  regions.  An  iterative  procedure  for  satisfying  the 
constraints  is  developed  which  can  be  applied  as  post-processing  to  standard  MMSE  beamformers.  The 
algorithm  is  used  to  synthesize  a  nearly  uniform  sidelobe  level  quiescent  pattern  for  the  UESA,  and  to  con¬ 
trol  sidelobe  levels  for  the  same  array  in  an  adaptive  manner.  Performance  results  using  data  provided  by 
Lincoln  Lab  show  that  under  low  sample  support  conditions,  sidelobes  can  be  effectively  suppressed  while 
maintaining  high  signal-to-interference  plus  noise  ratio,  and  deep  nulls  on  clutter  and  interferes. 
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1  Introduction 


Space-Time  Adaptive  Processing  (STAP)  used  in  airborne  radar  systems  combines  signals  from  N  antenna 
array  elements  and  M  pulses  to  adaptively  suppress  clutter  and  jamming  in  both  the  space  (angle)  and 
time  (Doppler  frequency)  dimension.  The  foundation  of  most  STAP  techniques  is  the  Minimum  Variance 
Distortionless  Response  (MVDR)  processor  [1],  The  standard  MVDR  processor  weights  are  designed  to 
minimize  the  processor  output  power  subject  to  a  linear  distortionless  constraint  in  the  angle-Doppler  steer¬ 
ing  direction.  The  same  weights  can  be  obtained,  to  within  a  scale  factor,  using  the  minimum  mean  square 
error  (MMSE)  criterion,  where  the  weights  are  designed  to  minimize  the  mean  square  error  (MSE)  between 
the  processor  output  and  a  reference  signal  [2]. 

Traditionally,  STAP  systems  have  used  a  rotating  linear  array  configuration,  however  a  fixed  circular 
ring  array  is  currently  under  development  under  the  UHF  Electronically  Scanned  Array  (UESA)  program 
sponsored  by  the  Office  of  Naval  Research  (ONR).  The  array  consists  of  54  directional  antenna  elements 
with  suppressed  backlobes  [3].  The  elements  can  be  modeled  as  having  a  cosine  pattern  as  shown  in  Figure 
1(a).  Only  a  portion  of  the  elements  will  be  used  for  each  steering  direction.  A  possible  configuration  is  to 
use  20  of  the  elements  at  a  time  to  transmit  and  receive  [3],  [4].  When  the  antenna  has  54  elements,  it  can 
be  scanned  mechanically  in  6.61°  increments  by  choosing  the  appropriate  20-element  sector,  and  scanned 
electronically  ±3.33°  with  the  chosen  sector  of  elements.  This  array  configuration  is  illustrated  in  Figure 
1(b).  The  resulting  array  spatial  beampattem  is  shown  in  Figure  2(a)  and  the  space-time  beampattem  is 
shown  in  Figure  2(b). 

The  circular  array  configuration  has  the  potential  to  provide  continuous  360°  availability,  however  it  has 
some  potentially  negative  impacts  for  STAP  algorithms  [4].  First,  the  clutter  rank  is  increased  in  a  manner 
similar  to  the  increase  from  misalignment  with  the  velocity  vector  in  linear  arrays.  Second,  the  clutter  locus 
varies  with  range.  This  decreases  the  number  of  range  gates  that  can  be  averaged  to  reliably  estimate  the 
clutter  covariance  matrix.  In  radar  systems,  this  results  in  large  sidelobes  in  the  MVDR/MMSE  beamformer 
which  can  lead  to  increased  false  alarms  from  clutter  and  unexpected  interferers.  Pattern  control  can  be 
improved  with  steering  vector  tapering  using  a  quiescent  pattern  with  desirable  main  beam  and  sidelobe 
characteristics  [5],  linear  main  beam  constraints  [2],  white  noise  gain  constraints  [6],  quadratic  quiescent 
pattern  constraints  [7]-[ll],  and  reduced  rank  subspace  techniques  [12]-[14]. 

In  this  project,  we  have  expanded  on  several  of  these  techniques  and  developed  a  general  framework 
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for  adaptive  and  non-adaptive  beampattem  synthesis  for  non-linear  arrays  based  on  minimum  mean-square 
error  (MMSE)  beamforming  with  quadratic  beampattem  constraints  (QPC).  In  this  technique,  main  beam 
and  sidelobe  pattern  control  is  achieved  by  imposing  a  set  of  inequality  constraints  on  the  weighted  mean- 
square  error  between  the  adaptive  pattern  and  a  desired  beampattem  over  a  set  of  angle-Doppler  regions. 
By  proper  choice  of  the  number  of  constraints,  the  angular  regions  to  which  they  apply,  and  the  desired 
beampattems  in  those  regions,  we  can  trade  off  the  level  of  pattern  control  with  algorithmic  complexity. 
At  one  extreme,  we  can  derive  low-complexity  beamforming  techniques  based  on  one  or  two  constraints 
similar  to  the  adaptive  pattern  control  methods  in  [7]-[10].  At  the  other  extreme,  we  can  achieve  tight  pattern 
control  using  many  constraints,  in  a  manner  similar  to  the  technique  in  [11].  A  more  moderate  approach 
using  several  constraints  is  shown  to  maintain  a  high  SINR  and  good  pattern  behavior  with  reasonable 
complexity.  The  algorithm  uses  an  iterative  procedure  for  satisfying  the  constraints  which  can  be  applied 
as  post-processing  to  standard  STAP  processors.  The  technique  can  also  be  used  for  non-adaptive  pattern 
synthesis.  It  generalizes  the  techniques  in  [15]  and  [11]  for  developing  low  sidelobe  quiescent  patterns 
for  arbitrary  arrays,  and  can  be  used  for  developing  the  desired  patterns  used  in  the  adaptive  method.  The 
technique  has  been  applied  to  the  UESA  for  both  adaptive  and  non-adaptive  pattern  synthesis.  Spatial 
beamforming  results  with  the  UESA  circular  array  were  published  in  [16]  and  CSTAP  results  will  appear  in 
[17]  and  [18]. 

This  report  is  organized  as  follows.  In  Section  2,  the  MMSE-QPC  problem  formulation  and  optimum 
solution  are  developed,  as  well  as  the  iterative  implementation.  In  section  3,  circular  array  STAP  results 
using  data  provided  by  MIT  Lincoln  Lab  [3]  are  presented.  A  summary  and  areas  for  future  research  are 
provided  in  Section  4. 

2  MMSE  Beamforming  with  Quadratic  Pattern  Constraints 

2.1  Problem  Formulation 

We  assume  a  STAP  model  with  N  antenna  elements  and  M  pulses.  Let  v(0,  </>,  u)  denote  NM  x  1  space- 
time  array  response  vector  to  a  signal  arriving  with  elevation  angle  9,  azimuth  angle  <f>,  and  normalized 
Doppler  frequency  u.  We  partition  azimuth  angle-Doppler  space  into  r  sectors,  fii, . . . ,  Qr,  as  shown  in 
Figure  3.  In  this  illustration,  the  elevation  angle  space  has  only  one  partition  and  the  sectors  are  cubes, 
however  more  general  partitions  of  azimuth  angle,  elevation  angle,  and  Doppler  space  can  be  used.  Let 
Bd,i  (9,  <f>,  u)  =  w %.v(9,  fa  u)  be  a  desired  beampattem  in  the  region  f 1,  ,  and  w d,i  be  the  corresponding 
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Figure  3:  Partition  of  Angle-Doppler  space 


weight  vector.  The  MSE  between  the  beampattem  generated  by  the  adaptive  weight  vector  w  and  the 
desired  beampattem  over  the  region  is  given  by 


€?  =  f  |whv(0,<?!>,w)  - 
The  error  can  be  written  compactly  as 


(1) 


e?  =  (w  -  Wj/Q.-fw  -  wd,i) 


(2) 


where 

Q i=  f  v(6,<t>)v(0,<j>)HdCl.  (3) 

Thus  the  pattern  error  is  a  quadratic  function  of  the  adaptive  weight  vector. 

We  design  adaptive  weights  according  to  the  standard  MMSE  beamforming  criterion,  while  limiting 
the  deviations  from  the  desired  pattern  using  quadratic  pattern  constraints.  The  MMSE-QPC  optimization 
problem  is 

min  w^Rxw  —  ■wHKxs  —  R^w  +  a\  (4) 

st.  (w  -  vrdti)HQi(w  -  w dii)  <  Li  i  =  1, . . . ,  r 

where  Rx  is  the  data  covariance  matrix,  RXi  is  the  cross-correlation  vector  between  the  data  and  reference 
signal,  and  a]  is  the  reference  signal  power.  When  the  reference  signal  is  the  desired  signal,  the  cross¬ 
correlation  vector  has  the  form 

(5) 
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The  solution  to  the  MMSE-QPC  optimization  problem  has  the  form 

w  =  ^Rx  +  ^2  {v2* v*  +  IZ  AtQtwdt,  j  .  (6) 

In  the  resulting  beamformer,  a  weighted  sum  of  ‘loading’  matrices  Q =  1, . . . ,  r  are  added  to  the  sample 
covariance  matrix,  and  a  weighted  sum  of  desired  weight  vector  terms  QjW^j, i  =  1, . . r  is  added  to 
the  steering  vector  vs.  The  loading  terms  balance  the  adaptive  pattern  with  the  desired  pattern.  There  are 
a  set  of  optimum  loading  levels  A,-,  i  =  1, . . . ,  r  which  satisfy  the  constraints,  however  there  is  no  closed 
form  solution  for  the  loading  levels,  even  when  r  =  1.  It  can  be  shown  that  the  mean-square  pattern  error 
decreases  with  increasing  A,-,  but  at  the  expense  of  decreased  interference  suppression.  The  loading  levels 
must  be  chosen  judiciously  to  achieve  pattern  control  while  maintaining  high  signal-to-interference-plus- 
noise  ratio  (SINR). 

When  the  array  is  a  uniformly  spaced  linear  array,  and  there  is  one  sector  which  covers  the  entire 
angular  space,  and  the  technique  is  the  same  as  diagonal  loading.  When  one  or  two  quadratic  constraints 
are  imposed  over  the  main  beam  and/or  sidelobe  regions,  beamformers  similar  to  those  developed  in  [7]- 
[10]  are  obtained.  These  are  relatively  low  complexity  techniques  that  have  been  shown  to  be  effective 
for  pattern  control  when  the  loading  levels  are  set  appropriately.  In  practice,  the  optimum  loading  level  is 
scenario  dependent  and  it  may  be  difficult  to  set  properly  to  achieve  the  required  pattern  control. 

When  there  is  a  single  main  beam  constraint  and  many  constraints  over  a  dense  grid  of  discrete  points 
in  the  sidelobe  region,  the  technique  in  [11]  is  obtained.  This  is  a  high  complexity  technique  that  achieves 
tight  pattern  control.  By  constraining  the  beampattem  at  a  set  of  fixed  points,  logical  constraint  levels  can 
be  easily  specified.  If  the  desired  pattern  in  the  sidelobe  region  is  set  to  zero,  the  pattern  error  is  just  the 
beampattem  value,  and  the  constraint  level  can  be  set  to  the  to  the  desired  sidelobe  level.  Constraints  may 
also  be  imposed  in  the  mainbeam  region  to  maintain  the  mainbeam  shape,  although  setting  the  constraint 
level  isn’t  so  straightforward. 

We  take  a  moderate  approach  using  several  constraints  to  obtain  good  pattern  behavior  with  reasonable 
complexity.  If  the  sectors  are  chosen  small  enough,  we  can  still  specify  the  desired  pattern  to  be  zero  in  the 
sidelobe  region,  and  enforce  maximum  sidelobe  levels  by  setting  the  constraint  value  to  the  desired  sidelobe 
level  times  the  volume  of  the  sector.  As  the  size  of  the  sectors  becomes  larger,  the  pattern  control  is  less 
exact,  however  the  complexity  is  reduced.  We  have  found  that  we  can  define  a  reasonable  number  of  sectors 
that  allow  us  to  specify  the  constraints  to  maintain  the  sidelobe  control  we  need,  at  a  significantly  reduced 


complexity  compared  to  the  technique  in  [1 1]. 

2.2  Iterative  Weight  Computation 

In  [11],  an  iterative  procedure  is  used  to  adjust  the  loading  levels  to  achieve  sidelobe  level  constraints.  A 
similar  procedure  can  be  used  to  obtain  the  appropriate  loading  levels  here.  The  algorithm  can  be  initialized 
with  the  standard  MMSE  weights 

w<°)  =  <ts2R^vsi  (7) 

or  with  tapered  adaptive  weights  obtained  by  replacing  a2vs  with  a  tapered  steering  vector  wt, 

w<°>  =  R^w*.  (8) 


At  each  iteration,  let  the  weight  error  in  each  sector  be 

„w  =  „<»-■>  -  wAj 


(9) 


The  pattern  errors  are  computed  and  checked  against  the  constraints.  If  a  constraint  is  exceeded,  i.e.  if 
w)^  Q,w^  >  Li,  the  loading  for  that  sector,  \\p\  is  increased  by  A-p^  to  A-p'  =  Ajp_1^  +  a\p\ 
otherwise  it  remains  at  its  current  level.  The  covariance  matrix  and  weights  are  updated  by 


Rxp) 

=  Rx  +  EAiP)Q. 

1=1 

(10) 

Wt(p) 

=  wt  +  EASP)Q«Wd,i 

1  =  1 

(ID 

w(p) 

(12) 

This  is  a  computationally  expensive  procedure  because  the  covariance  matrix  is  inverted  at  each  iteration. 
However,  if  the  incremental  loading  levels  a|p^  are  small,  the  update  can  be  accomplished  without  re¬ 
inverting  the  matrix  by  using  a  first  order  Taylor  series  approximation  to  the  weight  vector  wOO.  Let  S(p) 
denote  the  covariance  matrix  inverse  at  each  iteration, 

=  (r£>)  _1  =  (n(rl)  +  E  A1P)Q« 

\  i=l 

S(p-1))"1+EaJp)Q8)  .  (13) 

t=l  / 

The  steering  vector  at  each  iteration  can  be  written  as 

w<»>  =  w«'-»  +  £A!p)qiWAi.  (14) 

1=1 
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Using  (13)  and  (14)  in  (12),  the  adaptive  weight  vector  is  given  by 
wW  =  S^wjp) 

=  (»ri)+E4?w) 

The  first  order  Taylor  series  approximation  to  w(p)  has  the  form 


(15) 


(16) 


where  A(p)  is  the  r  x  1  vector  of  loading  increments,  is  the  derivative  matrix  evaluated  at  A(p)  =  0, 

'dwW  dw<p) 


jy(p)  — 
uo  — 


.5a[p)  c>a<p)J 

and  is  the  weight  vector  evaluated  at  A^p)  =  0, 


A(p)=0 


7lP)=W<P)|  =W(P_1) 


The  ith  column  of  the  derivative  matrix  is 

dw<p> 


I  A(p)  =0 


-1 


-1 


i=l 


aA«»>  =  -((s(’",))"  +  |>ip,QiJ  +Ea!',)q, 

((s<'-'>)'1+E4p,q.) 


+ 


Therefore 


dw(p) 


dA 


(p) 


A(p)=0 


=  -S(p  ^Q,-  (s(p  1lw|p  ^-w d,ij 


=  -S(p-X)Q i  (w(p-1l  -  wdii)  . 
Combining  (18),  (20),  and  (9)  with  (16)  yields 

w(p)  =  w^P_1^  —  S(P-1)  ^2  A^Qt'W^. 

t=l 

The  matrix  inverse  is  updated  is 

g(p)  _  s(p-!)  _  s(p_1)  A^QjS^p_1). 


(17) 


(18) 


(19) 


(20) 


(21) 


(22) 


1  =  1 
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A  simple  way  to  achieve  fast  convergence  while  ensuring  that  the  small  update  assumption  is  valid  is  to 
let  A^  be  a  fraction  of  the  of  the  current  loading  value,  i.e.  A^  =  aA-p\  where  a  in  the  range  0.3  to  0.5 
seems  to  work  well.  This  requires  that  the  initial  loading  level  be  non-zero.  One  possibility  is  to  initialize 
all  of  the  loading  levels  to  some  small  value  A0. 

In  summary,  the  algorithm  is  initialized  by 


Ai°>  = 

s<°)  = 

w(°)  = 

and  the  weights  are  updated  by 

1.  for  i  =  1, . . r 


if  w^  Q,w^  >  L{  then  A-p^  =  a\\p^  else  A^  =  0  (27) 

Ajp)  =  aJp_1)  +  Ajp)  (28) 

2.  QW  =  ^Ajp)Qi  (29) 

1=1 

3.  qW  =  ^Ajp)QlWip)  (30) 

1=1 

4.  w<p>  =  w(p-x)  -  (31) 

5.  S(p)  =  S(p_1)  -  S(p_1)Q(p)S(p_1).  (32) 


A0,t  = 


-1 


^Rx  +  Ao^QtJ 

S(0)  (wt  +  A0  - 


(23) 

(24) 

(25) 


1  =  1 


/rs 


2.3  Non-Adaptive  Pattern  Synthesis 

The  technique  can  be  used  for  non-adaptive  pattern  synthesis  by  letting  Rx  =  I.  It  generalizes  the  tech¬ 
niques  in  [1 1]  and  [15]  for  developing  low  sidelobe  quiescent  patterns  for  arbitrary  arrays,  and  can  be  used 
for  developing  the  tapered  steering  vector  used  in  the  adaptive  method. 


3  Results 


In  the  MIT  Lincoln  Lab  data  set  [3],  there  are  A  =  20  elements  and  M  =  18  pulses  with  a  300  Hz  pulse 
repetition  frequency.  First,  the  MMSE-QPC  technique  was  used  to  synthesize  a  -35  dB  uniform  sidelobe 
level  quiescent  pattern  steered  to  (f>  —  0°  and  u  =  0  Hz  for  a  range  of  50  km,  which  corresponds  to  9  — 
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Figure  4:  Partition  of  Angle-Doppler  space  for  UESA 

—  10.5°.  Angle-Doppler  space  was  partitioned  into  one  elevation  angle  sector  6  €  (—11°,— 2°),  11  azimuth 
angle  sectors  <fi  €  (-12°,  12°),  ±(12°,  30°),  ±(30°,  60°),  ±(60°,  100°),  ±(100°,  140°),  ±(140°,  180°), 
and  5  Doppler  sectors  u  e  (-30, 30),  ±(30, 90),  ±(90, 150)  Hz  for  a  total  of  1  x  11  x  5  =  55  sectors,  as 
shown  in  Figure  4. 

The  desired  pattern  was  set  to  zero  outside  of  the  mainlobe  region,  and  the  constraint  levels  were  set 
to  -35  dB  times  the  volume  of  the  sector.  No  constraint  was  used  in  the  mainlobe  region.  The  loading 
levels  were  initially  set  to  A0  =  10  dB,  and  then  iteratively  increased  in  the  sectors  where  the  constraint 
was  not  met  using  a  =  1.5.  The  exact  weight  update  was  used  because  this  procedure  can  be  done  off¬ 
line.  The  algorithm  converged  in  five  iterations.  The  initial  conventional  space-time  beampattem,  and  final 
-35  dB  sidelobe  level  pattern  are  shown  in  Figure  5,  along  with  several  patterns  at  intermediate  stages  of  the 
iteration.  The  final  -35  dB  sidelobe  level  pattern  steered  to  <j>  =  0°,  u>  =  60  Hz  is  shown  in  Figure  6. 

Next,  a  scenario  with  two  30  dB  interference-to-noise  ratio  (INR)  jammers  at  60°  and  -20°,  in  addition 
to  clutter,  was  considered.  An  8  km  training  window  (200  snapshots)  was  used  to  estimate  the  covariance 
matrix.  The  standard  MMSE  processor  weights  were  computed  by  adding  -30  dB  diagonal  loading  to  allow 
the  covariance  matrix  to  be  inverted.  The  resulting  space-time  beampattem,  and  beampattem  cuts  are  shown 
in  Figures  7(a)  and  7(b).  The  beamformer  has  put  nulls  on  the  clutter  ridge  and  the  two  jammers,  however 
the  sidelobes  are  quite  high. 

The  MMSE-QPC  adaptive  beamformer  was  used  to  reduce  the  sidelobes.  The  quiescent  weights  de- 
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Azimuth  Angle  (deg.) 


Azimuth  Angle  (deg.) 


(a)  Conventional  beampattem 


(b)  Initial  Loaded  Beampattem 


-90  0  90  180  -180  -90  0  90 

Azimuth  Angle  (deg.)  Azimuth  Angle  (deg.) 


(c)  Second  Iteration 


(d)  Fifth  Iteration 


Figure  5:  Quiescent  Pattern  Synthesis 


Doppler  Frequency  (Hz) 


-180  -90  0  90  180 

Azimuth  Angle  (deg.) 


(a)  Space-Time  Quiescent  Beampattem 


Beampattem  cut  at  60  Hz 


Beampattem  cut  at  =  0  deg. 


(b)  Quiescent  beampattem  cuts 

Figure  6:  Low  Sidelobe  Level  Quiescent  Beampattem  Steered  to  <f>  =  0°,  u>  =  60 Hz 
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rived  in  the  first  example  was  used  as  the  tapered  steering  vector.  The  initial  loading  levels  were  set  to 
Ao  =  0.0013,  and  then  iteratively  increased  using  a  =  0.8.  In  six  iterations,  the  MMSE-QPC  is  able  to 
reduce  the  sidelobes  below  the  -35  dB  level  while  maintaining  a  well  behaved  main-beam,  and  deep  clutter 
and  jammer  nulls.  The  final  beampattem  is  shown  in  Figures  8(a)  and  8(b). 

4  Summary  and  Future  Research 

The  MMSE-QPC  technique  provides  a  general  framework  for  quiescent  and  adaptive  space-time  beampat¬ 
tem  synthesis  which  provides  both  main  beam  and  sidelobe  control  with  reasonable  computational  com¬ 
plexity.  Beampattem  control  is  achieved  by  imposing  a  set  of  inequality  constraints  on  the  weighted  mean- 
square  error  between  the  adaptive  pattern  and  a  desired  beampattem  over  a  set  of  angle-Doppler  regions. 
An  iterative  procedure  for  satisfying  the  constraints  was  developed  which  can  be  applied  as  post-processing 
to  standard  MMSE  beamformers,  and  guidelines  for  setting  constraint  levels  and  algorithm  parameters  to 
achieve  desired  performance  were  discussed. 

The  algorithm  was  used  to  synthesize  a  nearly  uniform  sidelobe  level  quiescent  space-time  beampattem 
for  the  UESA,  and  to  control  sidelobe  levels  for  the  same  array  in  an  adaptive  manner.  Performance  results 
using  data  provided  by  Lincoln  Lab  showed  that  under  low  sample  support  conditions,  sidelobes  could 
be  effectively  suppressed  while  maintaining  high  signal-to-interference  plus  noise  ratio,  and  deep  nulls  on 
clutter  and  interferes. 

This  formulation  can  be  extended  in  a  straightforward  manner  to  the  MVDR  and  more  general  linearly 
constrained  minimum  variance  (LCMV)  beamforming  criteria.  It  can  also  be  used  in  other  fully  and  partially 
adaptive  STAP  architectures.  Topics  to  be  pursued  under  the  second  phase  of  the  project  include  applying 
the  quadratic  pattern  constraint  methodology  to  LCMV  beamforming;  generalizing  the  technique  for  fully 
and  partially  adaptive  STAP  architectures  in  element  and  beam  space,  and  with  pre-  and  post-Doppler  pro¬ 
cessing;  applying  the  technique  to  the  3D  STAP  model  for  hot  clutter  mitigation;  and  evaluating  subsequent 
adaptive  detection  performance. 
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Doppler  Frequency  (Hz) 


Azimuth  Angle  (deg.) 


(a)  Initial  adaptive  space-time  beampattem 


Be  am  pattern  cut  at  60  Hz 


Azimuth  Angle  (deg.) 


MVDR-GSC,  AG  =  16.8427  dB 


Doppler  Frequency  (Hz) 


(b)  Initial  adaptive  beampattem  cuts 

Figure  7:  Initial  Adaptive  Beampattem  Steered  to  <f>  =  0°,  u>  =  60  Hz. 
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Doppler  Frequency  (Hz) 


(a)  Adaptive  space-time  beampattem  after  sixth  iteration 


Be  am  pattern  cut  at  60  Hz 


(b)  Sixth  Iteration  beampattem  cuts 


Figure  8:  Final  Adaptive  Beampattem  Steered  to  <f>  =  0°,  u>  =  60  Hz. 
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